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Abstract 



We report the first fully general numerical calculation of the neutron and 
electron dipole moments, including the seven significant phases. We find that 
there are major regions in the parameter space where none of the phases 
are required to be small, contrary to the conventional wisdom. The electric 
dipole moments (EDM's) do provide useful constraints, allowing other regions 
of parameter space to be carved away. We keep all superpartner masses light 
so agreement with experimental limits arises purely from interesting relations 
among soft breaking parameters. 
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I. INTRODUCTION 



The general parametric structure of the Minimal Supersymmetric Standard Model 
(MSSM) includes a large number of CP-violating phases. Their presence has largely been 
ignored in phenomenological analyses because of severe constraints imposed on individual 
phases by the experimental upper limits for electron and neutron electric dipole moments if 
other phases are set to zero. These limits would generally constrain the phases considered 
individually to be less than 10~ 2 unless the mass parameters are pushed beyond the TeV 
scale UJ. Recently, it has been emphasized ||, however, that cancellations between different 
terms contributing to the dipole moments can allow for values of the phases very differ- 
ent from zero even when the superpartner masses are relatively light. Since this can have 
important consequences for extraction of the parameters in the SUSY Lagrangian from 
experimental data, for calculation of dark matter densities and scattering cross sections, for 
baryogenesis, for Higgs boson limits, and more, it is rather important to study the problem of 
constraints on the complex phases without making any unnecessary simplifying assumptions 
based on theoretical prejudice. To put it differently, the phases may or may not actually 
be small. We must find out from data, without making assumptions that lead to excluding 
regions of parameter space where parameters are large. If the phases are large, they affect 
many CP-conserving quantities throughout particle physics, so it is even more important 
to proceed carefully. The phases can only be large if certain approximate relations among 
soft breaking parameters hold; these relations will be checked in future experiments. The 
relevant relations are not fine tuned, but are quite reasonable, with various soft-breaking 
parameters related in size and sign to one another. 

Some important results have already been reported in literature. Nath and Ibrahim 
have presented @ some of the formulas needed in the analysis and calculated the EDM's 
in the framework of minimal supergravity model. Olive et al. JIJ have analyzed the case 
of two phases and confirmed the analysis of Nath and Ibrahim; they also applied it to the 
calculation of neutralino relic density and detection rate. Some similar results were reported 
in Ref. || and phenomenological consequences of CP-violating phases in the MSSM were 
studied in 0. 

In this paper we want to address this issue in its entirety in order to establish a connection 
between the usual parameters of the MSSM Lagrangian and ranges of the phases allowed 
by experimental data on the electric dipoles. We work in the framework of the simplest 
possible model neglecting the flavor mixing but avoiding any assumptions about unification 
of the soft breaking parameters. We use light superpartners masses, so apart from relations 
among soft breaking parameters the resulting EDM's would be very large. 

II. PHASE STRUCTURE OF THE FULL MSSM 

We define the MSSM to be the supersymmetric theory with the same particles as the 
Standard Model (SM) plus their superpartners, the SM gauge group, two Higgs doublets, 
and conserved R-parity. 

The MSSM Lagrangian |],||] depends on a total of 126 parameters and it includes three 
well known sources of CP-violating phases. The first is related to the two Higgs doublets 
present in the model since both the fi parameter in the superpotential and the soft breaking 
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parameter b can be complex and their phases are denoted and (fb respectively. Three 
more phases, <pi, ip 2 , and tp 3 , enter through the complex masses of the gauginos associated 
with the standard gauge groups. Finally, most of the phases originate in the flavor sector of 
the Lagrangian, either in the scalar soft mass matrices m;t _ ^ T _ or the trilinear matrices 

° ° ' Q,u,d,L,e 

a u ,d,e- The mass matrices are hermitian so only off-diagonal terms can be complex but the 
trilinear matrices are general 3x3 matrices allowing for the the diagonal entries to also be 
complex. 

The impact of the phases associated with the off-diagonal terms on experimental observ- 
ables is suppressed by the same mechanism which is required to suppress the existence of 
large flavor changing neutral effects and for the purposes of this study all these phases can be 
neglected; if some of them matter it will only strengthen our results. We assume that all the 
scalar soft mass matrices and trilinear parameters are flavor diagonal and that the complex 
trilinear terms are proportional to the corresponding Yukawa couplings a/ = yfAfe ipA f 1 
where ip^f are the relevant phases Q. 

It is important to realize that not all of the listed phases are physical. Above the elec- 
troweak symmetry breaking scale, the Lagrangian possesses two partial U(l) symmetries 
which can be promoted to full symmetries by treating the dimensionful parameters as spu- 
rions charged under those symmetries M. Under an R-symmetry the Grassmann variable 
9 (8) is charged +1 (—1) and therefore this symmetry distinguishes between component 
fields of the superfields. If the charge of a chiral superfield is r, its scalar component field <fi 
transforms under R-symmetry with charge r, the fermionic field if) has charge r — 1 and the 
auxilary F scalar field possesses charge r — 2. In order to preserve the R-invariance of the 
superpotential it is convenient to choose r = 1 for the matter superfields and r = for the 
Higgs superfields. The advantage of this choice is also clear from the fact that R-symmetry 
defined in this way is not broken in the process of electroweak symmetry breaking. The 
vector superfields are not charged under R symmetry, and so only the gaugino component 
field A (A) obtains a charge of +1 (—1). It is clear that to preserve R symmetry as a full 
symmetry of the superpotential and also of the soft SUSY breaking terms in the Lagrangian, 
lii Aj and the gaugino masses Mj have to be charged under the R-symmetry. 

The whole MSSM Lagrangian with the exception of the //-term in the superpotential 
and the b soft breaking term is also invariant under a Peccei-Quinn (PQ) symmetry. This 
symmetry transforms the Higgs fields with charge —2 and the matter fields Q, u, d, L, e with 
charge +1. Again, if \i and b are treated as spurions full symmetry is restored above the 
electroweak scale. Below this scale, the PQ symmetry is broken as the Higgs fields acquire 
vacuum expectation values. Physical observables can only depend on such combinations of 
parameters which are invariant under all symmetries of the Lagrangian. For the unbroken 
theory we have two symmetries and therefore two conditions, allowing us to eliminate two 
phases. When electroweak symmetry is broken we are left with only one unbroken symmetry, 
but the phase of b, which is related to the phase of the Higgs vacuum expectation values, 
can be absorbed into the physical Higgs fields by appropriate redefinition. It is therefore 
natural to take b (and the VEV's) to be real and set one more phase to be zero. We 



1 In further text we take (unless explicitly stated otherwise) all the dimensionful parameters to be 
real and positive; their phases are always written explicitly. 
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prefer in this paper to take ip2 = thus explicitly violating reparametrization invariance but 
one has to keep in mind that all other parameter choices are related to our choice by an 
R-transformation. A fully reparametrization invariant approach to CP-violation in SUSY 
theories will be discussed elsewhere || . Our numerical results do not depend in any way on 
this simplification. 

Taking into account our parametrization choice, the final set of set of phases considered 
in the discussion of the electron and neutron electric dipole moments includes three phases 
appearing in the chargino-neutralino-gluino sector, namely (fi, ip^ and (p^, and four phases 
(fA u , (pA d i y>At an d <fiA e corresponding to the trilinear soft breaking parameters relevant in the 
dipole moment calculation as discussed in the following section. As we will see below, even 
though a u , ad, a e are proportional to small Yukawa couplings, their phases enter because 
contributions to the EDM's require a chirality flip leading to dipole moments' proportionality 
to the relevant mass. 



III. ELECTRIC DIPOLE MOMENT CALCULATION 

The electric dipole interaction of a spin-1/2 particle / with an electromagnetic field is 
described by an effective Lagrangian 

Cedm = ~d f fa^ l5 fF^. (3.1) 

In theories with CP- violating interactions, the electric dipole df receives contributions from 
loop diagrams. The best way to account for such contributions is to use the effective theory 
approach in which the heavy particles are decoupled at some large scale Q and the full theory 
is matched with an effective theory including a full set of CP-violating operators If 



we restrict ourselves to dimension 5 and 6 operators the effective Lagrangian takes the form 

£*ff = Y,d{Q)Oi{Q) (3.2) 
i=i 

where the Ci(Q) are Wilson coefficients evaluated at scale Q, and the Oi are the three 
considered operators 

i _ 

Ox = -~f a^ls fFp,, (3.3) 

2 = ~fWsT a fG^, (3.4) 

3 = -±f abc G« p G b iGl a e^. (3.5) 

(3.6) 

It is obvious that all three operators contribute when the external fermionic particles are 
quarks, while in the case of the electron Cf and Cf are identically zero. 

Supersymmetric models contribute to the Wilson coefficients at the one loop level and 
they include several types of graphs as shown in Fig. 1. Chargino, neutralino and gluino 
loops where the second particle in the loop is a scalar superpartner, either a slepton or a 
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squark, contribute to G\ and G% coefficients depending on whether a photon or a gluon is 
radiated. The contributions can be calculated at the electroweak scale since a typical SUSY 
scale in most models is of the same order of magnitude. For the gluino loop contribution to 
the quark EDM the matching gives 

Cr\Q) = - j^EMAr^i^), (3.7) 

cr\Q) = ^EMAT s )^c(^). ( 3 - 8 ) 

The neutralino and chargino loops contribute both to the electron and quark electric dipole 
moments and one finds 



C(k ~ N{Q) = ^^ Qf ^ lm{A %~ N) ^ B ^ (3 ' 9) 
477 sin Uw i=1 J=1 m 4 rrij 



n 2 6 4 - Tn ~ m> 2 ~ 

cr N (Q) = ff^EEMAr")-* ) (3.io) 

1U/I i=l j=l ' i "H 



and 



ct- c (Q) = — . 2 E E )— t[Q/-b(— t) + (<?/ - tf/)^(-4)], (3-11) 

cr c (Q) = -f^EEMAr )^^)- (3-12) 

1077 2 ~~J~~L m « m i 

In equations 3.7-12, are the masses of the corresponding scalar particle running in the 
loop and A, B and C are the loop functions obtained by integrating out the heavy particles 
in the loop. These functions, together with the vertex A functions calculated in our phase 
parametrization, can be found in the Appendix. Qf denotes the electric charge of the 
external fermion and Q'* is the charge of internal sfermion when different from Qf. 

The gluonic operator O3 obtains a contribution from the top-stop loop with a gluino 
exchange as shown in Fig. 1 and one has 

f oq n - m 2 — m~ ml m 2 m 2 

C((Q) = -3a s m t (^) 2 lm(A u r 9 ) h 5 t2 H(-±, -|, ^|), (3.13) 
477 m5 m~ m~ mi 

where the loop function H can be found in the Appendix. 

The Wilson coefficients then have to be evolved from the decoupling scale Q down below 
the chirality breaking scale A x using the renormalization group equations (RGE's) in order 
to account for resummation of the logarithmic corrections. So far, only QCD corrections for 



quark operator Wilson coefficients have been estimated |L2| for RGE evolution down from 
the electroweak scale to A x giving 

C?(A x )=r li C?(Q), (3.14) 
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where rji ~ 1.53 and 772 — — 3-4. All other corrections are neglected in our calculation. 
At the low scale, the CP-violating operators Oi have to be projected on the electric dipole 
operator to evaluate their contribution to the numerical value of the electric dipole. This 
is a complicated task since the chirality breaking scale A x = 1.18 GeV is very close to the 
QCD scale and perturbative methods are not reliable in this region. The best thing one can 
do at present is to use naive dimensional analysis |13| which yields 



d f = C{(A X ) + ^Qf(Ax) + ^C/(Ax)- (3-15) 

Finally, since the neutron is a composite particle, one has to use a phenomenological neutron 
model to calculate the neutron EDM from the moments of the constituting quarks. From 
the simple SU (6) quark model one obtains 

d n = ^(4d d -d u ), (3.16) 

where d d and d u are the EDM's of the down and up quark respectively. 

One of the important features of the contributions to the EDM is the fact that the 
effective Lagrangian in Eqn. 3.1 requires different chirality of the the initial and final particle. 
In the supersymmetric diagrams this can happen in two ways — either the exchanged squark 
or slepton change chirality via L-R mixing terms in the sfermion mass 2 matrices and couple 
to the gaugino component of the intermediate spin- 1/2 particle, or the L and R sfermions 
preserve their chirality and couple to the higgsino components of charginos or neutralinos. 
As a result, all contributions are directly proportional to the mass of the external particle 
since both the L-R mixing sfermion mass term and the higgsino-fermion-sfermion coupling 
are proportional to the relevant Yukawa coupling. Another consequence of the chirality flip 
is the explicit proportionality of the contributions to the mass of the intermediate spin-1/2 
particle. 



IV. CONSTRAINTS ON THE PHASES 



As already mentioned, we present a numerical treatment of the electric dipole moment 
calculation, with the main emphasis on the cancellations between various contributions to 
the Wilson coefficients. This allows large values of the phases to give contributions consistent 
with the experimental bounds on the values of the electric dipole moment of both the electron 
and the neutron. Current experimental limits for the neutron require that 
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IcLl < 1.1 x 10 



-25 



ecm 



(4.1) 



and for the electron JT5[ 



\<L\ < 4.3 x 10 



-27 



ecm, 



(4.2) 



at 95% confidence level. 

We start our analysis by choosing a simple set of the MSSM parameters which leads 
to a fairly light spectrum of the superpartners while still keeping the general set of the 
seven relevant CP- violating phases, tpi, tp 3 , ip^, <Pa u , ^PA d , <PA t an d ^A e , which we consider. 
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Since our results do not assume heavy spectrum suppression of the CP-violation effects, 
they are fairly general in the sense that increasing the masses of the superpartners can 
only broaden the effect of cancellation between different contributions to the electric dipole 
moment. The resulting ranges of the phases for different spectra will differ quantitatively 
from our examples but the general observation that the phases indeed do not have to be 
small will still remain valid. To simplify the set of parameters we assume that the squark 
and slepton soft masses and the trilinear soft parameters are flavor diagonal. We also take 
the diagonal entries of these matrices to be universal for all three generations and neglect 
any splitting between up-type and down-type right-handed squark masses. Similarly, for 
the sneutrinos we consider a single universal soft mass for all three flavors. The trilinear 
parameter A is assumed to be not only flavor universal but also the same for both sleptons 
and squarks. As a result, we are left with the following set of soft parameters in the scalar 
flavor sector — mc,, m^, m~ lR rriq L , m,g R and A. We do not assume any relation between 
the gaugino masses other than taking M 1 < M 2 < M 3 . Unless stated otherwise, all of our 
calculations consistently employ a common set of parameters shown in the following table: 



Standard set of parameters 
(values at EW scale) 


Mi = 75 GeV 


m 9 = 185 GeV 


M 2 = 85 GeV 


m h = 195 GeV 


M 3 = 250 GeV 


m in = 225 GeV 


H = 450 GeV 


nig L = 340 GeV 


A = 250 GeV 


vriq R = 360 GeV 


tan/? = 1.2 


m A = 300 GeV 



We have varied them sufficiently to show that our results are qualitatively unchanged 
for significant regions of these parameters. 

A. Electron EDM 

Let us now concentrate on discussing the cancellation mechanism in the two cases of 
EDM calculation. The electron EDM limits are more constraining than the neutron EDM 
limits and are also simpler to study since there are only two large contributions. The electron 
EDM calculation involves only the chargino and neutralino contribution to the C\ Wilson 
coefficient corresponding to the electric dipole moment operator and only three phases enter 
the calculation, namely ip^, tpi and ipA e - Since the neutralinos are mixtures of both U(l) 
and SU(2) gauginos and both neutral higgsinos, the neutralino contribution includes both 
types of chirality flipping processes — from the gaugino exchange with L-R slepton mixing 
as well as from the process with gaugino-higgsino mixing and requiring no chirality flip in 
the slepton sector. On the other hand, the chargino exchange can only proceed through 
the latter channel since the chargino SU (2) gaugino component only couples to left-handed 
fields. In order for a cancellation between these two contributions to occur, certain conditions 
have to be met. 

The first condition requires that the two contributions have opposite sign over at least 
a subset of the phase parameter space. In fact, this requirement is automatically satisfied 
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for contributions coming from the gaugino-higgsino mixing diagrams. These contributions, 
involving both charginos and neutralinos, depend on tp^ and have opposite sign over the 
whole range of ip^ due to the fact that the \x parameter enters the neutralino and chargino 
mass matrices with opposite phase. This "fortunate" feature can be traced back to the 
antisymmetry of the SU{2) metric e appearing in the superpotential. 

The neutralino contribution of this type could in principle also depend on ipi, which would 
upset the exact anticoincidence of the signs. In practice, the dominant part of the contribu- 
tion comes from the SU (2) gaugino-higgsino mixing and the effect of ipi is constrained to a 
shift in the phase of the neutralino contribution from the gaugino-higgsino interaction. The 
gaugino-gaugino diagram for the neutralino contribution, unlike the gaugino-higgsino dia- 
grams, involves L-R mixing in the selectron sector. The imaginary part of the the relevant 
phase dependent term is 

lm(—^N*]) ~ 2 m e[^ e sin((^ 1 + (p Ae ) + ^tan/3sm(<^ - c^)], (4.3) 

m LL m g 

where m? LR and rn? LL are elements of the selectron mass matrix. In order for this expression 
to have the oposite sign to the chargino contribution, which is negative for < ip^ < tx and 
positive for n < ip^ < 2n, various possibilities occur depending on the relative sizes of A e 
and \x tan (3. For example, in the two limiting cases where |/Lt| tan /3 >> \A e \ and \A e \ » 

tan/3, we get in terms of the phases that modulo 2n we have to impose ip^ — ipi ~ — (p^ 
and <pi + pa e ~ — W respectively. In general this contribution is opposite in sign to the 
chargino one over a significant part of parameter space. 

The second condition of cancellation requires the chargino and neutralino contributions to 
be of the same magnitude. In the chargino sector, the gaugino-higgsino mixing involves only 
one type of gaugino while in the neutralino sector there are two gaugino states and therefore 
the elements of the neutralino diagonalizing matrix N generally yield smaller imaginary 
parts than the elements of the chargino matrices U, V. Moreover, the chargino contribution 
is enhanced due to larger values of the A(x) loop function as compared to the B(x) function 
in the neutralino expression. This comes from the fact that the photon in the chargino 
loop diagram is emitted from the fermionic leg of the corresponding diagram as opposed to 
the scalar leg in the neutralino diagram. Both these effects decrease the relative magnitude 
of the neutralino contribution compared to the chargino contribution. On the other hand, 
in the neutralino case the gaugino-gaugino contribution can balance some of the difference 
between the two contributions. For that to happen, it is important that the relative size 
of the chargino-higgsino contribution decreases and the relative size of the gaugino-gaugino 
contribution increases, which can be achieved by increasing /i. That brings us back to the 
first condition of opposite sign which can be satisfied for tan/5 >> \A e \ if p x ~ n. 

As a result, for suitable combinations of the dimensionful parameters an almost exact 
cancellation can occur for the whole range of ip^ as exemplified in Fig. 2 a. In this plot we 
chose fi = 700 GeV in addition to our standard set of parameters and ipi was set to be 
equal to ir. The values of tpA e were varied randomly leading to the result that the values of 
the neutralino contribution and of the total dipole moment form bands of non-zero width, 
while the chargino contribution is independent of (p Ae - It is clear from the plot, however, that 
virtually all values of ip^ would be allowed for this particular set of parameters depending only 
on a suitable choice of the ip Ae value range. This is also significant because ip Ae is otherwise 
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irrelevant not only in the neutron EDM calculation but also in most other phenomenological 
considerations. Later figures show effects of varying \x. Note that without cancellations one 
would have to have each contribution reduced by ~ 1CT 2 , i.e. each phase would have to be 
< 1CT 2 , as in the usual result. 

Thus we see that the strong constraint from the electron EDM limit is naturally satisfied 
over a significant part of the parameter space, though not all of it. The neutralino and 
chargino contributions can automatically have opposite sign and the same magnitude for 
most of the ip^ range when the mass parameters are in certain ratios depending on the 
other phases. While the cancellations do require related magnitudes of some parameters, 
and could thus be interpreted as a fine-tuning, we think that the required mass relations are 
typically the kind of relations that might arise in a theory of the soft-breaking parameters, 
and are likely to be a clue to the form of the theory. The resulting relations are predictions 
that can be tested in other experiments. 

To put it in another way, there are two ways to satisfy the electron EDM constraints. 
One possibilty is that the phases are small or zero as a result of some presently unknown 
mechanism. Alternatively, the phases could be large, and the masses could have certain 
approximate ranges of reasonable values. The relevant signs would automatically give the 
needed cancellation, which need not have happened. The two alternatives lead to very 
different predictions for many other observables. The naturalness of the cancellation that 
occurs leads us to consider the solution with large phases seriously enough to convince us 
to analyze the full parameter space and to study the resulting predictions, which we will 
report on later. 



B. Neutron EDM 

Next we turn to the neutron EDM, where cancellations are easier to obtain. First of 
all, all three operators in Eqn. 3.2 receive contributions from the MSSM one loop diagrams 
involving quarks as incoming and outgoing particles. The gluino-squark diagram projects 
on both Oi and O2 operators, and the contribution of the relevant Wilson coefficients Cf k ~ g 
and C| fc ~ 9 to the EDM is numerically comparable. The contribution of C| fe_ff is seemingly 
suppressed by the factor of ^ in Eqn. 3.15 compared to Cl k ~ g , but that is compensated by 
enhancements from the factor of ^f, and mainly from the loop function C(x) in the matching 
conditions 3.7 and 3.8. This is again a consequence of the fact that the gluino leg in the 
diagram can emit gluons but not photons. The chargino loop C{ k contribution is typically 
of the same order as the gluino loop contributions while C^ - *" contributes negligibly since 
in this case the — enhancement alone does not overcome the suppression from j-. Both 

neutralino contributions from Cf k ~ N and C| fe_7V can be safely neglected in the neutron dipole 
analysis. Reasons similar to those for the electron case lead to the suppression of Cf k ~ N 
compared to Cf k ~ c ', but in the quark case this effect is more pronounced since the squarks 
are typically heavier than the sleptons and they have fractional charges. Correspondingly, 
the contribution from Cl k ~ N is even smaller than that from Cf k ~ N . This effectively reduces 
the number of phases by eliminating (pi as one of the parameters numerically relevant in the 
neutron EDM calculation. 

As in the electron case, it is necessary that the chargino contribution be opposite in 
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sign to the sum of the other three contributions for the cancellation to occur. The gluino 
contribution exhibits the same behavior as the gaugino part of the neutralino contribution 
in the electron case and Eqn. 4.3 transforms into 

Im^G* 2 ) ~ -\m q \A q sin(y? 3 + ^0 + /*/(/?) Mv, ~ Vs)], (4.4) 

m LL m e 

where A q = A u , Ad and f((3) = cot f3, tan/? for up and down type quarks respectively. The 
contribution from the pure gluonic operator , on the other hand, depends only on p 3 , ip^ 
and A t as far as phases are concerned and the role of (p^ and pA t is again determined by 
the relative size of fi cot j3 and A t . This implies that p> 3 and tp^ are the crucial phases in 
the EDM calculation. In order to demonstrate the cancellation on a practical example, in 
Fig. 2b we set p 3 equal to ir and take all three trilinear parameter phases to be consistent 
with zero. All these choicess are enforced within a small variation around the central value 
leading to a non-zero width of the gluino and pure gluonic contribution. In addition we 
choose \i = 300 GeV so that it is comparable in magnitude to A q = 250 GeV and the off 
diagonal squark mixing terms get a comparable contribution from both terms in Eqn. 4.4. 
The resulting sum total of the neutron EDM is consistent with zero over a wide range of 
(p^. As (p 3 , <pa u , ^A d and (pA t are varied this situation will persist for large but correlated 
ranges of these phases. The variation with fi is shown in later figures. 

C. Numerical results 

The effects of the cancellation mechanism on the ranges of phases allowed by the EDM 
experimental limits can be explored by varying all phases randomly for a given set of mass 
parameters and plotting the allowed points projected on planes in the phase parameter 
space. In Fig. 3a and 36 we show the allowed regions for the standard parameter set with 
H = 450 GeV in the — <p\ and p^ — <p 3 plane respectively. The filled black circles signify 
the points allowed by electron EDM constraints and the open circles stand for those allowed 
by the neutron EDM limits. The ipi dependence has little significance as far as the neutron 
constraints are concerned, while in the electron case tpi has to be correlated with ip^ in order 
to satisfy the limits. Only a selected band of the values of ip^ is allowed by the electron 
constraints and the neutron constraint imposes a correlation between the values of p 3 and ip^ 
within this band. Still, when these conditions are satisfied, values of very different from 
or 7r are allowed. All values of p\ and ip 3 can occur while the EDM limits are respected. 

In Fig. 4 we display the same results as in Fig. 3 but we take /x = 60 GeV. The range of 
(Pfi is constrained by both electron and neutron limits in this case, and the interval allowed 
by both is significantly narrower than in the previous case. Nevertheless, all values of <pi 
and (p 3 are permitted again. 

It is important to see how the range of allowed values of ip^ depends on ji because 
iffj, plays a crucial role in the electron as well as in the neutron EDM calculation. Fig. 5 a 
displays this range for both calculations with the standard parameter set and varying /1. The 
overall trend shows that for larger values of \i it is easier to satisfy the EDM limits. The b 
frame shows the effects of A e variation on the electron EDM constraints when ji = 450 GeV. 
Similarly, in Fig. 6 we examined the dependence of the allowed ip^ range on an overall 
scaling parameter x which rescales all the dimensionful parameters in the standard set and 
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\i = 450 GeV according to the formula M' = xM. It is interesting to note that in order to 
allow the full range of ip^ one has to go to fairly large parameters x > 4 while the same 
effect can be obtained by raising /i to be larger than 450 GeV. 

Finally, in Fig. 7a and b we plot the lightest neutralino mass vs. ip^ and ipi respectively 
for the standard parameter set and \x varied from 50 GeV to 800 GeV. The neutralino masses 
can vary quite dramatically in the allowed regions and this fact substantially affects phe- 
nomenological observables at colliders and cosmological implications of the supersymmetric 
model. 

V. CONCLUSION 

We have shown that the role of the cancellation mechanism in the calculation of the 
electron and neutron electric dipole moments within the general framework of the MSSM 
including a non-restricted set of CP-violating phases has crucial consequences for the range 
of individual phases. Even with a light sparticle spectrum, phases can have values very 
different from zero and still satisfy experimental bounds on the values of the electron and 
neutron EDM's. 

A trivial but possible way to avoid constraints from the dipole moment measurements is 
the traditional one that all supersymmetric phases are equal to zero or unnaturally small. 
This would require the existence of some presently unknown mechanism which would ensure 
that there is negligible CP violation in the SUSY breaking sector of the MSSM Lagrangian. 
On the other hand, we have found that the phases may be large while certain approximate 
relations hold among the mass parameters and phases, resulting in cancellations in the 
calculation of the electron and neutron EDM. These relations could in principle also come 
from a theory of SUSY breaking predicting the exact form of the soft SUSY breaking sector 
in the Lagrangian. We have shown in this paper that the latter possibility is legitimate and 
the ultimate decision between the two alternatives should be made based on experimental 
measurements. 

We have presented a study of the constraints imposed on the phases by electron and 
neutron EDM data for some particular values of soft parameters with relatively light spec- 
tra. The results exhibit general features typical for similar choices and they show that all 
considered phases can have non-zero values, ip^ is severely constrained while other phases 
can have any value as long as certain correlations with ip^ are respected. The constraints 
on the phases relax as heavier spectra or large values of /i are considered. 

The fact that phases can be non-vanishing is very important if one considers the gen- 
eral correspondence between the parameters in the supersymmetric Lagrangian and various 
observables which will possibly be measured at future collider experiments. For example, 
without a determination of the phases it is not possible to measure the value of tan (3. It 
is also important to realize that the presence of phases has a substantial impact on the 
neutralino relic density calculation and on the magnitude of the corresponding neutralino 
scattering cross section for dark matter detection. If progress in supersymmetric particle 
physics proceeds by the historical path, it will be essential to measure the phases to learn 
the form of the soft-breaking Lagrangian, and thereby be led to recognize the mechanism of 
supersymmetry breaking. 
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APPENDIX 



In this appendix we summarize all the calculational details necessary for evaluation of 
the contributions to electric dipole moments of elementary particles in the MSSM at the 
one-loop level. The effect of CP-violating phases enters through the particular vertex A 
functions characteristic for each type of contributing diagrams. These functions depend 
on the matrices diagonalizing the generally complex mass matrices of the participating 
super symmetric particles. 

In our parametrization, the gluino mass is complex and can be diagonalized by a single 
complex number G defined by 



G*M 3 e i ^G~ 1 = M 3 , 
resulting in G = e llfi3 ^ 2 ''. Simlarly, the chargino mass matrix 

' M 2 



M 



a 



y/2M w sin (3 

itp/j, 



^ V2M W cos (3 ixe 
is diagonalized by two generally complex unitary matrices U and V so that 



WMcV- 1 = M d c ag . 



(A.l) 



(A.2) 



(A.3) 



The neutralino mass matrix contains two phases, (fi and <£> M , in our parametrization 
/ 

M N = 



Mie w -M z sin 9 W cos (3 M z sin 9 W sin (3 \ 

M 2 M z cos 9 W cos f3 — M z cosOw sin/3 

M z sin 9 W cos (3 M z cosO w cos(3 -/xe^" 

V M z sin 9 W sin (3 -M z cos9 w sm(3 -fie^ 



(A.4) 



and the diagonalization matrix satisfies 

N~ U M N N = M 



N 



(A.5) 



Finally, the scalar superpartners of the three families of fermions in the Standard model 
obtain masses through a general mass matrix 



m Q,Q,L + mu.d^m^ e + D L 1 a l ude v uAd - fxe l ^v diU , a l 
a Ujd , e Vu,d,d - fie-^Vd^ul m 2 - A - + m], d e m Uidie + D R 1 



t 



M 



u,d,e 



(A.6) 



where D L = M Z (T 3 — Q sin 9^) cos 2/3 and D R = M Z Q sin 9^ cos 2/3, and v u , v d are the 
VEV's of the two neutral Higgs fields coupling to the up-type and down-type particles 
respectively. The parameters in bold print are 3x3 matrices, generally complex as discussed 
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in the main text. The mass matrix in Eqn. A. 6 can be diagonalized by a pair of 3 x 6 matrices 
relating the interaction and mass eigenstates 

//• = r , /,; ; /f" ? ' (A.7) 

f? = r'jv, /f" ?; (A.8) 

for each type of fermion and all families i = 1, 2, 3. Our notation distinguishes between the 
three types of sfermions, u d and e, and individual flavor states are numbered according to 
the family number, so, for example, the u% field corresponds to the left-handed top squark 
field. 

The gluino vertex function reflects the fact that the gluino is a pure gaugino and the only 
possible way to produce a chirality changing effective vertex is to make use of L-R squark 
mixing and get 

i — 1 (q)ki L (q)ki U l A - y J 

with no summation implied over i and q = u,d. The neutralino vertex function can be 
obtained in a similar way giving 

A f k -N> = { ^2 tan e wQN * jT f q)kl - \ f N* hj rf q)H }x 

{-V2[tan MQ " T 3 )N% + T 8 JVyrg w + N* hj T? q ] ki } (A.10) 

Wllere A - = v^C'sin/T Ad ' e = V2M^cos^ aIld h = 3 fOT h = aIld h = 4 fOT f = U - lt 18 

obvious from the structure of the function that the neutralino effective vertex includes both 
gaugino and higgsino interactions. Finally, the chargino vertex function for individual types 
of particles takes the form 

A«* J = ^uV*2Tfd)ki(UjlT(d)ki ~ ^dU* 2 ^fd)ki) (A.ii) 
A»* ' — ^dU* 2 ^f u )ki(yjl^fu)ki ~ ^dV* 2 ^fu)ki) (A- 12) 

Ar dj = AeC/; 2 y; x . (A.13) 

In order to make this paper self contained, we also list the necessary loop functions 
coming from integrating out the supersymmetric particles in the one loop diagrams in the 
case of the electric and chromoelectric dipole operators |TJ| 

1 , 2ln(x) 
2(1 -xf K 1-x 

C(I)= l^ (Iflx _26 + ^=W-if^) (A.16) 
o(l — x) z 1 — x 1 — X 

and from the two loop calculation in the case of the purely gluonic operator [T7] 

1 f l f l f l N X N 2 
H{zx,z 2 , z 3 ) = 77 / dx du dyx(l-x)u , (A. 17) 

I Jo Jo Jo 



A ( x ) = oTi 1^(3 ~x+ -^r) (A.14) 
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where 

Ni = u(l — x) + z 3 x(l — x){l — u) — 2ux[ziy + ^2(1 — y)] 

AT 2 = (1 - a;) 2 (l - uf + u 2 - \x 2 (l - uf 

9 

D = u(l — x) + z 3 x(l — x)(l - u) + ux[z x y + z 2 (l - y)]. 
The integrals in the above definition of H can be simplified and evaluated numerically. 
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FIGURES 



FIG. 1. One loop Feynman diagrams contributing to the calculation of the electric dipole 
moments in the MSSM. The gluon and photon line can originate on any internal leg carrying 
corresponding charge. 

FIG. 2. Illustration of the cancellation mechanism in the EDM calculation. See the discussion 
in the text. Frame a includes the contributions to the electron dipole moment arising from neu- 
tralino and chargino loops contributions to the C\ Wilson coefficient for varying ip^, ipi ~ and 
values of tpA e sampled randomly. A standard set of parameters (see text) with \i = 700 GeV was 
used. Frame b shows the neutron EDM contribution from the gluino loop graph projection into C\ 
and C2 {§1 and g~i, from the chargino loop contribution to C\ (C±) and from the gluino-top-stop 
graph contributing through the purely gluonic operator Wilson coefficient C3 (G). In this case, 
the standard set of parameters is adopted with \i = 300 GeV, ip% ~ -it and pA q ~ for q = u,d,t. 
In both cases the natural cancellations can give a total of order the experimental limits for most 
or all of ip^ . If the cancellation effects were not included one would conclude that all phases would 
have to be of order 10 -2 to not exceed the experimental limits. 

FIG. 3. Plots of regions allowed by the electron (filled circles) and neutron (open circles) EDM 
limits in the ip^ — <pi plane (frame a) and the (p^ — ps plane (frame b). A value of fi = 450 GeV 
was chosen together with the standard parameter set and all phases were sampled randomly. 

FIG. 4. Same as Fig. 3, but for /j, = 60 GeV and the standard set of parameters. Again, all 
phases were varied randomly. 

FIG. 5. Frame a shows variation of the (p^ allowed region with n for the standard set of 
parameters and other phases sampled randomly. The values of A = A e = A u = A^ = At were 
also varied from —500 GeV to 500 GeV. Open (full) circles denote points allowed by the neutron 
(electron) EDM limit. Frame b demonstrates variation of the p^ range allowed by the electron 
EDM limits with the values of A for ^ = 450 GeV. 

FIG. 6. We plot the points allowed by the electron EDM limits for parameter sets with all the 
mass parameters scaled by x with respect to the standard set. All phases are sampled randomly. 

FIG. 7. Plots of the lightest neutralino masses allowed by the neutron (open circles) and 
electron (filled circles) EDM limits vs. ip^ in frame a and ipi in frame b. In addition to all phases, 
the values of \x were also varied from 50 GeV to 800 GeV and all other parameters were standard. 
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Fig. 1 
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Fig. 3 
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Fig. 6 
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